—— a 


МАСА TN 3150 9976 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 3150 


METHOD FOR RAPID DETERMINATION OF PRESSURE CHANGE 
FOR ONE-DIMENSIONAL FLOW WITH HEAT TRANSFER, 
FRICTION, ROTATION, AND AREA CHANGE 
By James E. Hubbartt, Henry O. Slone, and Vernon L. Arne 


Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


Washington 
June 1954 


ШІ 


WN ‘УЯ лнунап нда. 


і КШ... 


CU-1 


TECH LIBRARY KAFB, 


ЇШШЇШЇЇЇ 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 5150 


METHOD FOR RAPID DETERMINATION OF PRESSURE CHANGE FOR ONE-DIMENSIONAL 
FLOW WITH HEAT TRANSFER, FRICTION, ROTATION, AND AREA CHANGE 


By James E. Hubbartt, Henry 0. Slone, and Vernon L. Arne 


SUMMARY 


An approximate method for rapid determination of the pressure 
change for subsonic flow of в compressible fluid under the simultaneous 
action of heat transfer, friction, rotation, and area change is de- 
veloped. Іп the development of this method, the momentum equation was 
approximated and rearranged for a convenient solution employing charts. 
This report presents both the analysis involved in simplifying the mo- 
mentum equation and the charts necessary for obtaining particular solu- 
tions. The charts provide a step-by-step solution which converges to 
an exact solution as the number of steps is increased, An illustrative 
example and comparison with more rigorous numerical solutions with con- 
ditions typical for air-cooled turbine blades are included. These com- 
parisons show that the solution converges rapidly to provide good 
accuracy. 


INTRODUCTION 


The effective design of ducts to accommodate air flow requires раг- 
ticular solutions of the momentum equation for determining the pressure 
changes encountered. Simplified methods which give an accuracy suffi- 
cient for engineering purposes are in demand. An approximate method 
for determining such particular solutions has been developed for the 
one-dimensional flow of a compressible fluid under the influence of heat 
transfer, friction, rotation, and area change. 


A number of studies of one-dimensional flow of a compressible fluid 
have been published. Reference 1 presents a rather complete treatment 
including a form of the momentum equation with heat transfer, friction, 
and area change which is suitable for numerical integration. As an aid 
for expediting the numerical integration, coefficients of the differ- 
entials in the momentum equation of reference 1 have been tabulated as 
"influence" coefficients in reference 2. The analyses of references 1 
and 2 have been extended in reference 3 to include rotational forces. 
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Reference 5 also provides a form of the momentun equation suitable for 
numerical integration and tabulates "influence" coefficients for ex- 
pediting the calculations. However, in all cases, the numerical inte- 
grations as employed in references 1 to 3 are tedious and time consum- 
ing. Particular solutions suitable for rapid determination of pressure 
changes for one-dimensional flow with both arbitrary and specific heat- 
input distributions have been presented in references 4 to 7. However, 
these solutions include only heat transfer and friction. In addition, 
reference 8 developed an approximate solution for constant flow area in- 
cluding rotation as well as heat transfer and friction. The approxima- 
tions of reference 8 were based upon experimental data for cooled tur- 
bine blades tested under static conditions. Although the time required 
for a particular solution of the momentum equation and subsequent eval- 
uation of the pressure change was appreciably reduced when the approxi- 
mate solution of reference 8 was utilized, the time required for a large 
number of particular solutions, as frequently required for design pro- 
cedures, remains very lengthy. In addition, reference 8 does not account 
for flow-area changes frequently encountered in actual practice. 


Since a method for the rapid determination of the pressure change 
through ducts with heat transfer, friction, rotation, and area change 
is in demand, the momentum equation has been approximated and rearranged 
for a convenient solution employing charts. The purpose of the present 
report is to present the analysis involved in simplifying the momentum 
equation so that a graphical solution is possible and to present the 
charts necessary for obtaining particular solutions for subsonic flow. 
Heat transfer, friction, rotation, and area changes are all included in 
the solutions. The charts provide a step-by-step solution which con- 
verges to an exact solution as the number of steps is increased and is 
limited only by the graphical accuracy. Comparisons are also made with 
the solutions determined from the more exact numerical integration of 
reference 3 for typical examples of air-cooled turbine blades. 


ANALYSIS 


The general momentum equation for one-dimensional flow is approxi- 
mated and rearranged to form the basis for the construction of charts 
suitable for the determination of the Mach number change with flow under 
the simultaneous action of heat transfer, friction, rotation, and area 
change. A relation between the Mach numbers and pressures is finally 
determined such that charts may be employed to evaluate the pressure 
changes іп а duct once the Mach number changes are known. The energy 
changes may be arbitrarily specified in terms of the total-temperature 


changes and may be evaluated by в procedure such as that given in ref- 
erence 5. 
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Evaluation of Mach Numbers 


The differential form of the general momentum equation for one- 
dimensional flow is (ref. 3) 


ar 


ар - Fog dx + pV ат «200 (1) 


where х increases in the flow direction (all symbols are defined in 
appendix A). The differential drag force can be expressed by the 
equation 


2 2 2 
аға = f BV as = ғ РЎ 7 ах = 26 У д ах (2) 


2 2 Т, 2 


and the body forces per pound of weight flow, as resulting from rota- 
tion about a central axis, can be expressed by the equation 


2 
т 
eor i 


If the differential momentum change is replaced by the equivalent ех- 
pression 


pv av = alov?) - v a(pv) (4) 


and equations (2) to (4) are substituted into equation (1), the momentum 
equation, as integrated between stations x and х+Ах, becomes 


Pax i (oV?) лу (OV) Ax 
dP + alov?) - v а(рт) - 
Px (ov2).. (оу) 


х 


x-+Ax хълх 2 
werp ах + <= om ax = 0 (5) 


Оһ 


x x 


The first two terms of equation (5) сап be integrated directly.  How- 
ever, the last three integrals of equation (5) cannot be evaluated 
directly and therefore are approximated in order to obtain forms which 
can be conveniently solved. The third integral, which exists only in 
the case of area change, and the fifth integral, which expresses the 
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drag forces, шау be conveniently approximated by using the integrand ав 
the average of the values at the end points. That is, 


(oV). Ax (V. +y ) 
V alov) = XE Гру) ады - (ОТ) | (6) 
(ғу), 
and 
x+Ax 


4Е 


2 4f 2 2 
£V" ах м бк NEST. (e) + х e Й (7) 
2 Dy xx х+Ах Раух 


Іп order to facilitate a convenient solution, the fourth integral, 
which expresses the rotational forces in equation (5), will be approxi- 
mated by either 


2 
Dh 


x 


x+Ax . 
|, аўтр ax = ағу, Ах (8а) 


ог 


х+Ах | 
| авто dx s Or AP Axa (8b) 
x 


depending upon whether the solution is initiated with the conditions at 
the station x or the station x+\x. Since the integrand of either 
equation (8a) or (8b) is approximated by the.value at a single end point, 
the approximation is in general less accurate than those given by equa- 
tions (6) and (7). Equations (6) to (8) become: identities in the limit 
as Ax approaches zero, 


Because the density and velocity variations from stations x to 
x+Ax are approximated by equations (6), (7), and (8), the energy equa- 
tion must be employed only to determine the total energy change. 


When the approximations given by equations (6), (7), and (8a) are 
substituted into equation (5) and the first two terms of equation (5) 
are integrated, the approximate momentum equation becomes 
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Чу + Талк 
Руди = Px + (0V*) ang = (7x - Am Гот) ал = (ox - 
ағ 2 4f 2 
ی‎ [eem (2) ЛЕ 


However, if equation (8b) is used instead of equation (8a), the approxi- 
mate momentum equation becomes 


у +V 
Бал - Рк + (o72) - (072), - Ea EAE [(от) зд - (OV) | - 


4f 2 4f 2 | 
2 Ах х+Ах (ex x (x ) 
D TAP جي‎ + | m i + گے‎ | L— =0 9b ) 
x х 2 БЕЗ 2 x+Ax Dn x 2 ( 


If the variables of equations. (9a) and (Sb) are changed by employing 
the continuity equation, the relation between the total and static tem- 
peratures, and the equation of state, the two forms of the approximate 
momentum equation can be reduced to the useful relations (appendix B) 


ғыда Aca) = Ме (Е) С + TuS sto | 

10а) 

ала | 

Fite) = PRE (e [rns - ханз ЕТ) 
10b) 

where | 
(MA) = ey? fa - ач) (ax) 
cop = 002220) i (a2) 


gRM 
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ағ,Ак АМ 
A, = Bax - 422 - 1) (13) 


The total-temperature ratio in equations (10) may be evaluated from 
a solution of the energy equation such as presented in reference 3. 
Evaluation of Pressures 


By applying the continuity equation and the equation of state, the 
weight flow can be expressed as 


_P _ PVreR 
W = пи VA = nyt МА (15) 
or 
PA 42/2 
"= Це м) (16) 


Since the relation between the total and static temperature as given 
by the energy equation is 


Seite . (17) 


| 


equation (16) can be rewritten in the form 


1/2 
XT. TE (2 + ze)" (18) 


Using P'/P = (т'/туї/(т-1) to express equation (18) in terms of the 
total pressure gives 


Eu) 
WT ви (1 (riae) rA (19) 
P'A R 2 
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Equations (18) and (19) are plotted in figure 1 for ү = 1.40 with 
wVT'/PA and w T'/P'A ав the ordinates and М as the abscissa. 
From the Mach number the value of either vat, /Р'А ог wa T /PA can 
be determined from figure 1. Then from the given conditions and area, 
Р or Р! can be computed. 


CALCULATION PROCEDURE 


Equations (10a) and (100) represent the final forms of the momentum 
equation to be used in solving for the Mach number at either station x 
or хаАх. If the conditions are known at station x, the Mach number at 
station x+tAx is determined from equation (10a); while if the conditions 
аге known at station х+Ах, the Mach number at station x is determined 
from equation (10b). Increasing values of X, corresponding to в posi- 
tive Ax, must be in the flow direction. A given duct can be divided 
into an arbitrary number of increments in Ax (not necessarily of the 
same length), and the Mach number change through that duct can then be 
computed by progressing stepwise from the end where all conditions are 
known. The size of each increment should depend upon the accuracy de- 
sired. Subsequent examples will indicate the convergence of the solu- 
tion as Ax is decreased. 


Both Е(М,А) and G(M) have been plotted in figure 2 for air with 
Y = 1.40 to expedite the solution. Since y has a secondary effect 
on the Mach number change and the variation of y encountered in air 
flow is small, the consideration of other values of у is unwarranted. 


The procedure for determining the Mach number change and thus the 
pressure change with air flowing in a duct is to use equation (10a) or 
(105) and figures 1 and 2. For example, the procedure for determining 
the Mach number at station x from the conditions at station x4^x is 
to first evaluate A yay (eq. (14)). From the value of the Mach number 


at station xx and A,iay, the values of F(MyjAx, Лулу) and 
с(М лук) are read from figure 2. Then, with the specified values of 


Or AXAX/'Pis x. and МЛ (А„/А алх), Р(М,,А,) is calculated 
from equation (10b). From F(M,,A,) and А, (ea. (15)), the Mach num- 


ber at station x is read from figure 2. Once the Mach numbers are 
known, the pressure change can be determined from figure 1. To illus- 
trate this procedure, a typical sample calculation for an air-cooled 
turbine blade including heat transfer, friction, rotation, and area 
change has been carried out in table I with the passage divided into 
two equal parts. Double columns (colum 9 and columns 20 to 27) are 
included in the sample calculation setup in table I во that the calcu- 
lations could proceed from either the x or xx station. For each 
of these double columns, the left column is applicable if the solution 
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is initiated at station х+Ах, while the right column is applicable if 
the solution is initiated at station x. The first 11 columns of table 
I give the conditions which must be prescribed. Although the friction . 
coefficient f is assumed to be the same at both the x and xx | 
stations, it may be allowed to vary if greater accuracy seems warranted. 
The remaining columns of table I give the computations required to de- . 
termine the Mach number at either station x or xx. The sample cal- 
culation shown is for the determination of the Mach number at station 

x from conditions at station хх, 


RESULTS AND DISCUSSION 


Іп order to evaluate the accuracy of the solution presented herein, 
calculations were made for conditions typical of those occurring for 
air-cooled turbine blades (1) by use of the charts and (2) by the method 
of numerical integration presented in reference 5. The results of these 
calculations employing equation (10b) with 1, 2, and 4 steps with incre- 
ments of equal length are tabulated in table IT and compared with the 
results obtained from the more detailed numerical methods of reference 
3. Three cases are illustrated in table II: (I) heat transfer and 
friction; (II) heat transfer, friction, and rotation; (III) heat transfer, 
friction, rotation, and area change. For the first two cases, the air 
temperature was assumed to vary linearly from the inlet to the exit. 
However, for the third case, the air temperature was varied exponentially 
to correspond with flow inside a duct heated by crossflow and with con- 
stant inside and outside surface heat-transfer coefficients. The total 
pressures used for evaluating the percentage differences between those 
obtained by the two solutions were computed from figure 1 and the known . 
conditions and areas. 


It is shown in table II that the solutions converge rapidly. In 
addition, the comparison between the two solutions shows that good agree- 
ment is obtained by using only two steps. Even for the third case, the 
computed total pressures agreed within 4 percent for only two steps. In 
this case the disagreement was essentially eliminated by using four steps. 
Solutions were also made with eight steps, but for the conditions used 
in the examples of table II, the improvement over the use of four steps 
was probably not within the accuracy of the charts. The approximations ` 
given by equations (6) to (8) are least accurate at the higher flow Mach 
numbers where slight changes in the total pressure cause large changes 
in the Mach number (fig. 1) and, therefore, large changes in the veloc- 
ity and state conditions. For this reason, the Mach numbers used in the 
examples of table II are intentionally high. It should be pointed out 
that the accuracy of the approximations employed in equations (10) are 
dependent upon Mach number changes, and not the passage length, when го- 
tational forces are neglected. Thus, the accuracies shown by the results 
of table ІІ are indicative of the accuracies which would be obtained for 
longer passages with similar Mach number changes due to heat transfer, 
friction, and area change alone. ” 
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In many cases the air temperature varies exponentially with x as 
in case III of table ІІ. Іп order to check the effect on the Mach num- 
ber change of using the more easily obtainable linear air temperature 
distribution instead of an exponential variation, the examples presented 
in table II were recomputed as follows: (1) Cases I and II, which were 
calculated for table II with a linear air temperature distribution, were 
determined for an exponential distribution, and (2) case III was com- 
puted for a linear variation instead of the exponential air temperature 
distribution used for table II. The results of these calculations showed 
that the differences between the two variations were negligible and prob- 
ably of the same order of magnitude as the accuracy of the charts. For 
this reason, it appears that although the air temperature might vary ex- 
ponentially with x, for convenience в. linear variation may be assumed 
for evaluating the temperatures needed for the incremental subdivisions. 


CONCLUDING REMARKS 


A method is developed for the rapid determination of the Mach num- 
ber change and pressure change for the flow of a compressible fluid in 
а duct with heat transfer, friction, rotation, and area change. The 
method uses a step-by-step integration which, when applied to specific 
examples typical for air-cooled turbine blades, converged rapidly and 
provided good accuracy with, at most, only four steps. In addition, for 
these specific examples, a negligible difference occurred between using 
а linear or an exponential temperature variation along the duct length. 
This result indicates that the more convenient linear temperature vari- 
ation may be used for evaluating the temperatures for the incremental 
subdivisions if the actual temperature varies exponentially, as fre- 
quently occurs in practice. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 14, 1954 
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APPENDIX А 


SYMBOLS 


The following symbols are used in this report: 


flow erea, sq ft 
duct length or span, ft 


hydraulic diameter, 4A/1, ft 
generalized body force, Ib/ib of coolant, w’r/g (eq. (5)) 


drag force, 1b 


2 
М\1 + 2 М 


friction coefficient 


_ 3/2 
(2 + Ll м) / 


it (eq. (12)) 
standard acceleration due to gravity, 52.2 ft/sec? 
wetted perimeter, ft 

Mach number relative to passage 

static pressure, 1ъ/за ft abs 

total pressure with ‘respect to passage, 1b/sq ft 
gas constant, 55.5 ft-1b/(1b)(9R) 

radius, ft 

surface area, sq ft 


Static temperature, OR 


total temperature with respect to passage, OR 
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ү velocity relative to passage, ft/sec 
” weight flow, 1b/ sec 
х distance measured from passage inlet in flow direction, ft 
Y ratio of specific heats, 1.40 
4f 
A. x _ Да - 1) (eq. (15)) 
Dy x Ax Ax 
AP xx Ay. 
x +AX 
Алу - R———— - )ر‎ - ) (еа. (14)) 
D 
h, x4+Ax Ах 
р mass density, slugs/cu ft 
Q аёулх 
ща 
eb angular velocity, radians/sec 
Subscripts: 
e exit of passage 
3 inlet of passage 


x distance measured from passage inlet in flow direction, ft 
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APPENDIX В 


DERIVATION ОҒ APPROXIMATE MOMENTUM EQUATION 


Upon collecting terms and employing the relations 


(eV). (=) (pV), (оу) алх 


Vx 
(oV). Ax (=) = (OVD Ax (oV ш > tev), 


and _ 


(от), (ак) = (өт), СР как (“ах _ (912) (OV) 


2 (COR 2 (РУ Ах 


equation (9a) becomes 


(eV?) as (oV), AE AA 

Presse = Pa = 2-2 ULL үтте 
(ov?) (БУ) ли 48, Ах 

e PM 


and equation (9) becomes 


GV e | LE. „Е хадхак 


2 
a = EE 2 Cm ۽‎ == "n - Or pax PxA = 0 


If the continuity equation is employed to obtain the relation 
= (VA), = (PVA) Ay 


or 


(оу), _ Ахлу 


ov х+Ах А, 


(20a) 


(20b) 


(21) 
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and if the simplifying notations 


4f Ay 
EE Consi) (15) 
4f AAE Ах+дх 
х+Ах “Бшм - збе ] (14) 
or Ax 
2-°5 
are used, equation (20a) becomes 
2 y2 
Ручак 7 Px ЕЕ (Ал - 4) , CT (Ay - 4) -% 0,50 
(22a) 
and equation (20b) becomes 
(ov?) ‚ (97 ) 
Poux T Рх 7 ЕА л, T 4 7 A, 7 43-8, Ax хх - 0 
(226) 


A more convenient form of equations (22a) апа (22b) сап be obtained if 
the variables are changed by employing the equation of state, the energy 
equation, and the continuity equation. From the equation of state, 

pV can be expressed as : 


ov? = В. 2. PM (23) 


In addition, the static pressure can be eliminated by using the con- 
tinuity equation and the equation of state as follows: 


w= БОЛУ = == АУ (24) 


or 


(25) 
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where the total- to static-temperature ratio can be eliminated by using 
the energy equation in the form 


І L-L y2 - 
т-ізіе-м | (17) 
to give the expression 
Дара Б 
м( + Et ë) 
2 
Finally, application of the continuity equation gives 
gpAV gA V 


tp = Ж E- Tl 
TU? < ي‎ ETE (28) 


or, by incorporating equation (17) in equation (28), 


(2 «та м) 
Тр = Ear 2. (29) 


If equations (23), (26), and (29) are substituted into equations (22а) 
and (22b) and wi/R/yg is factored out, the two forms of the approxi- 
mate momentum equation become 


a-j E CED S,‏ ف 
X. TR ш‏ - - ---- 172 
n‏ کک mage)‏ 
(30a)‏ 


and 
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аё " 
ES dog -TE суал E 
(30b) 


Upon simplifying, rearranging, and employing the definition of &, equa- 
tion (30a) becomes 


x afr 
F (Mex Қол) Ез (Sees) oes + т эб (102) 


and equation (30b) becomes 


т! А afr , Ax 
(Mg A = (еды hua) - тр eae) 
х \Ахых Ах 


(10b) 
where 
F(M,A) = i Е - тм (ко | (11) 
7. - X8 4 
м + rz м2) 
e0 - (a А ri а)“ 


(12) 
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TABLE I. - SAMPLE CALCULATION AND SETUP 


DC I ЕСТІСЕ ЕЕЕ 


Ax, | Ay, AxaAx? | Юу, Бай ар iw Толу x| 9, АА лу | Ах 
ft | aq ft q ft ft ft |rt [roues (3)/(4) A, Fag =, 26, 
1: e) C 


x [soa osos ann nn [oco то аю ESRI SEATE ETS 
Балана e aola [vn лв sas la ool a 


Mex | | Малы) | GCM) 
х105 x10? 


[eae | xem | eee ис pw зис 
E -0.5655 а ад | ша | овие [аво | | 


дау | 


тоати ) вм, A. Ем, А) + 
аё afr 
in = TR ек), 
(23)+(24) 


ш 
| 


Ак | Tx 
ТА NU 
(15)(14)2/2 


OSTE ML VOVN 


LT 
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TABLE ІІ. - COMPARISON OF RESULTS FROM NUMERICAL AND 
GRAPHICAL SOLUTIONS 
Inlet Mach 


number My 
from charts 


0,419 (0.406 |0. 41610. 418| 0.0310 
.505 . 
.141) .125 


"Heat transfer and friction. 
Heat transfer, friction, dnd rotation. 


“Heat transfer ; friction, rotation, and area change. 
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Mach number, M 


Figure 1. - Evaluation of welght-flow parameter from Mech number 


(ratio of specific heats, 1.40). (А 155- by 17-іп. print of 


this Pig. is attached.) 
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Figure 2. - Evaluation of F(M,A) and G(M) (ratio of specific 


heats, 1.40). 
attached.) 


(А 12- by 20-in. print of this fig. is 


рін,АЙ 
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13.0 
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(b) Mach 


Figure 2. - Continued. 
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attached.) 
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Mach number 


number range, 0.075 to 0.25. 


Evaluation of F(M,A) апа G(M) (ratio of 
(А 145- by 20-in. print of this fig. is 
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IN 


ANS 


on 
ЕҢ 


E зо 


+25 430 -35 А +45 
Mach number 


(с) Mach number range, 0.20 to 1.0. 


Figure 2. - Concluded. Evaluation of F(M,A) and G(M) (ratio of 


specific heats, 1.40). (A ui by 20-in. print of this fig. is 


attached.) 
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